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We have calculated numerically the gap exponent of the anisotropic Heisenberg model in the 
presence of the transverse magnetic field. We have implemented the modified Lanczos method to 
obtain the excited states of our model with the same accuracy of the ground state. The coefficient of 
the leading term in the perturbation expansion diverges in the thermodynamic limit (A'^ — > oo). We 
have obtained the relation between this divergence and the scaling behaviour of the energy gap. We 
have found that the opening of gap in the presence of transverse field scales with a critical exponent 
which depends on the anisotropy parameter (A). Our numerical results are in well agreement with 
the field theoretical approach in the whole range of the anisotropy parameter — 1 < A < 1. 
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I. INTRODUCTION 

The study of transverse magnetic field on the low di- 
mensional spin systems has been attracted much inter- 
est recently from experimental and theoretical point of 
view. The experimental observations^'^ on the quasi-one- 
dimcnsional spin-1/2 anti-fcrromagnet CS2C0CI4 is a re- 
alization of the effect of non-commuting field on the low 
energy behavior of a quantum model. This shows a quan- 
tum phase transition from the spin-flop phase (ordered 
antiferromagnetically in y-dircction) at the low magnetic 
filed to a paramagnet for high fields. Moreover a novel 
behavior has been observed in the specific heat close 
to the quantum critical point^. A connection between 
the ground state properties of the anisotropic Heisen- 
berg model (XXZ) in the transverse field to the reported 
quantum phase transition has been given by the quantum 
renormalization group approach^,. In addition, a recent 
mean field approach for the weakly coupled chains^ has 
given a very good agreement on the finite temperature 
phase diagram of the two dimensional model with the 
observed experimental data-'^. 

The spin-(s = ^) Hamiltonian of the XXZ model in a 
transverse field on a periodic chain of N-sites is 



N 



sfsy^,+Astst+,-hs^], (1) 



where J > is the exchange couplings in the XY easy 
plane, — 1 < A < 1 is the anisotropy in Z direction and h 
is proportional to the transverse field. The spin (s") on 
site n is represented by ^a^;a = x,y,z in terms of Pauli 
matrices. This model is a good candidate for explaining 
the low temperature behavior of CS2C0CI4. However, 
since the integrability of XXZ model will be lost in the 
presence of transverse filed, more intensive studies from 
the theoretical point of view is invoked. 

When h = 0, the XXZ model is known to be solvable 
and critical (gapless)^. The Ising regime is governed by 
A > 1 while for A < —1 it is in the ferromagnetic phase. 



Magnetic field in the anisotropy direction commutes with 
the Hamiltonian at /i = and extends the gapless region 
(quasi long range order) to a border where a transition to 
paramagnetic phase takes place. The model is still inte- 
grable and can be explained bya conformal field theory 
with central charge c = 1 (Ref.ja] and references therein). 

Adding a transverse field to the XXZ model breaks the 
U{1) symmetry of the Hamiltonian to a lower, Ising-like, 
which develops a gap. The ground state then has long 
range anti-ferromagnetic order (—1 < A < 1). However 
due to non-zero projection of the order parameter on the 
field axis it is a spin-flop Neel state. In fact at a spe- 
cial fleld {hci = a/2(1 -I- A)) the ground state is known 
exactly to be of the classical Neel type«£^. The gap van- 
ishes at the critical field he, where the transition to the 
paramagnetic phase occurs. Classical approach to this 
model reveals the mean field results ~ which is exact at 
s — > 00. The implementation of the quantum renormal- 
ization group*^ shows that the transition at he is in the 
universality class of the Ising model in Transverse Field 
(ITF). The phase diagram of XYZ model in transverse 
field has been also presented in Ref.Q- The scahng of 
gap, phase diagram and some of the low excited states 
at hci of the XXZ model in the transverse field has been 
studied in Ref. 10]. In this approach the scaling of gap is 
given by the scaling of operator s^ which is read from the 
asymptotic form of the correlation functiorkii {{sfs'f^^)). 
This correlation function contains two terms, an oscillat- 
ing and a non-oscillating one. Each part defines a specific 
scaling exponent for s^ which depends on the anisotropy 
parameter (A). 

Exact diagonalizationi?: and Density Matrix Renor- 
malization Group (DMRG)i^ results give us some knowl- 
edge on this model but not on the scaling of gap. A 
bosonization approach to this model in certain limits 
leads to a nontrivial fixed point and a gapless line which 
separates two gapped phases^^, moreover the connection 
to the axial next-nearest neighbor Ising model (ANNNI) 
has been addressed. The applicability of the mean-field 
approximation has been studied by comparing with the 



DMRG results of magnetization and structure factor-'^^. 
Recently the effect of longitudinal magnetic field on both 
the Ising model in Transverse Field (TF)i& and the XXZ 
model in TF has been discussedii. 

Here, we are going to present our numerical results on 
the low energy states of the XXZ model in the transverse 
field which has been obtained by the modified Lanczos 
method introduced in Sec. II. It is believed that in this 
approach the accuracy of the excited state energies is the 
same as the ground state energy. In Sec. Ill, we will ap- 
ply a scaling argument presented in Ref.[l3 to our model 
and examine it by the numerical results to obtain the gap 
exponent. We have then applied a perturbative approach 
in Sec. IV to study the divergent behavior of the coeffi- 
cient in the leading term of the perturbation expansion. 
The gap exponent has been obtained by the relation to 
the exponent of the diverging term. The gap exponent 
depends on the anisotropy parameter in agreement with 
the field theoretical results^*'. Finally, we will present the 
summary and discussion on our results. 



II. MODIFIED LANCZOS METHOD 

The theoretical investigation of numerous physical 
problems requires an appropriate handling of matri- 
ces of very large rank. Even if in many applications 
the matrix is sparse, the problem can not be solved 
by means of a direct diagonalization by standard rou- 
tines. The Lanczos method and the related recursion 
methodsi^iSSiSiiSS, possibly with appropriate implemen- 
tations, have emerged as one of the most important com- 
putational procedures, mainly when a few extreme eigen- 
values (largest or smallest) are desired. Grosso and Mar- 
tinelli have presented a relevant implementation of the 
Lanczos tri-diagonalization schemeS^, which allows to ob- 
tain a very fast convergence to any excited eigenvalue and 
eigenfunction of H, overcoming memory storage difficul- 
ties. To explain this method briefly, let us consider an 
operator H, with unknown eigenvalues Ei and eigenfunc- 
tions \tpi)- Any auxiliary operator A — f{H) commutes 
with H, and thus shares with it a complete set of eigen- 
functions corresponding to the eigenvalues Ai = f{Ei). 
In order to obtain the nearest excited state of H to any 
apriori chosen trial energy Et, we consider the auxiliary 
operator A in the form A ~ {H — Et)"^ . In a completely 
different context, this form is suggested by numerical 
analysis-^' to solve the schrodinger equation within any 
desired energy. Now, we are faced with the solution of 
the following eigenvalue equation 



A\^,) = {H-Etf I V»)=A.|V'^) 



(2) 



Our strategy to solve Eq.® is based on the Lanczos al- 
gorithm. We briefly summarize some basic features of 
the Lanczos procedure in its standard formulation. 
Let us denote with Lpi {i = 1,2, ..., A^) a complete set of 
bases functions, for the representation of the operator H 
(and hence of A). Starting from a seed state \uq), given 



by whatever chosen linear combination of the (pi, a set 
of orthonormal states |'Uo)l'Wi), ■•■, |wjv) is constructed via 
successive applications of the operator A as follows 



\Ui) = {A-ao)\uo), 
ao = {uo\A\uo). 



(3) 



In general 



\Un+l) = A\Un) -an\Un) -bn\Un-l), n > I. (4) 

The (non-normalized) state \Un+i) allows us to deter- 
mine the coefficients bn+i and a„+i of the (n-|-l)th iter- 
ation step, via the procedure 



bl+i 



Mi+i 



= {Un+l\Un+l) 

{Un+l\A\Un+l) 

{Un+l\Un+l) 



(5) 



where 60 = is the initial condition. After normalization 
of the state |u„+i) = j-^!— |C/„+i), the steps and Q 
are repeated with n replaced by n+1. In the new basis 
\un), the operator A is represented by a tridiagonal ma- 
trix Tjn, whose elements a„ and 5„ are explicitly known 
for (m < N) . The diagonalization of tridiagonal matrix 
Tjn gives the eigenvalues, Ai = {Ei — Et)"^. 

The transformation to the tridiagonal matrix is trun- 
cated at some stages because of the round-off error. How- 
ever, the ground state energy can be obtained up to some 
significant digits. The accuracy of the excited energies 
is lost in the usual Lanczos method {A = H) by the 
round-off error. The modified Lanczos method explained 
above allow us to get the higher energy levels with the 
same accuracy of the ground state energy. We can se- 
lect the tuning parameter Et in the range of accuracy of 
our method to get the excited energy levels. By choosing 
the appropriate Et , we got the energy gap of the model 
presented in Eq.QJ up to 8 digits. 



III. THE SCALING ARGUMENT AND GAP 
EXPONENT 

The XXZ model is integrable and its low-energy prop- 
erties are described by a free massless boson field the- 
ory. In the transverse magnetic field, Dimitricv et.al 
considerediS, the perturbed action for the model as 



S = Sa + h f dtdxS''ix,t), 



(6) 



where 5*0 is the Gaussian action of the XXZ model. The 
time dependent correlation function of the XXZ chain for 
A |< 1 has the following asymptotic formii 



(5^(x,t)S'^(0,0)) 



{-irA, 



(x2+w2t2)9/2 

A2 
' trj.2 _j_y2^2\{e/2+l/2e)' 



(7) 



where Ai and A2 are known constants^^, t 
imaginary time, v is the velocity and 



1 



arc cos A 



it is the 



(8) 



The exponent of energy gap has been estimated^" by us- 
ing the long-distance contribution of the oscillating part 
to the action via the scaling of S^ 



gap^ h'', 



M 



1 



1-9/2' 



(9) 



If the scaling of 5*^ is read from the non-oscillating part 
of {S^{x,t)S^ (0,0)) the contribution to the action gives 
the following scaling for the energy gap?^:^ 



gap' 



4-6-1/ 



(10) 



The smaller value of /x and i' defines the leading order 
of the dependence of gap on the transverse field, the gap 
exponent. Thus, /i is the gap exponent for —1 < A < 
and 1/ for < A < 1. 

We have implemented the modified Lanczos algorithm 
on the finite size chains (N=8,10,12,...,24) by using pe- 
riodic boundary conditions to calculate the energy gap. 
We have computed the energy gap for different values of 
— 1 < A < 1 and the chain lengths. The energy gap as a 
function of the chain length (N) and the transverse field 
{h) is defined 



gapiN,h) ^ E2iN,h) - EoiN,h), 



(11) 



where Eq is the ground state energy and E2 is the sec- 
ond excited state one. The first excited state crosses the 
ground state N/2 times for a finite chain and the last 
crossing occurs at the classical point^S, hd — ^2(1 4- A). 
These two states form a twofold degenerate ground state 
in the thermodynamic limit where Ei — Eq vanishes. 

In the case of h equals to zero, the spectrum of the XXZ 
model is gapless. The gap vanishes in the thermodynamic 
limit proportional to the inverse of chain length. 



N 



lim gap{N, h = 0)- 



B 



(12) 



The coefficient B is known exactly from the Bethe ansatz 
solutions^. We consider this equation as the initial con- 
dition for our procedurai^. Adding the transverse field to 
the Hamiltonian, a nonzero gap develops. The presence 
of gap can be characterized by the following expression, 



gap{N, h) 
9ap{N,0) 

in the combined limit 



TV- 



CX), 



l + /(a:), 



0, 



(13) 



(14) 



where x = Nh^ is fixed and f{x) is the scaling function. 
Thus, the gap at finite h can be defined 



gap{N, h) 



B 

N 



B 

N 



/(^) 



B 

N 



9{NW 



(15) 




Nh 




Nh 



FIG. 1: The product of energy gap and chain length [Ngap) 
versus Nh^ . (a) A = 0.5 (b) A = —0.5. A linear behavior 
is obtained by choosing e = 2.0 for all different chain lengths 
A^ — 18, 20, 22, 24. The solid lines are guides for eye. 



where g{N) is a function of only N . It is imposed that the 
function g{N) approaches a non-zero constant value in 
the thermodynamic limit: limAr^oo 9{N) — constant = 
Coo ■ The regime where we can observe the scaling of gap 
is in the thermodynamic limit {N — > 00) and very small 
value oih {h <^ 1). This means that the scaling behavior 
is observable at large x {x — Nh'^ ^1). The asymptotic 
behavior of g{N) ioi N —^ 00 defines that the large x 
behavior of f{x) must be proportional to x, 



I{^) 



9{N) 
B 



C 
lim /(.t) = -^x. (16) 



Thus, if we consider the asymptotic behavior of f{x) as 

f{x)^xt (17) 

the (/)-exponcnt must be equal to one (0=1). If we 
multiply both sides of Ea. (|15l) by N we get. 



Af- 



lim Ngap{N, h) = B 

+ aD(a;>l) 



Co 



(18) 



Ea. (|18|l shows that the large x behavior of Ngap{N, h) is 
hnear in x where the scahng exponent of the energy gap 
is e. 

We have plotted the values of Ngap{N, h) versus 
Nh" for A = 0.5,-0.5 in Fig.|Tl). The reported re- 
sults have been computed on a chain of length N = 
18, 20, 22, 24 with periodic boundary conditions. Accord- 
ing to Ea. (|18|l . it can be seen from our numerical results 
presented in Fig.l^a) and Fig.(^b) that the linear bea- 
havior is very well satisfied by e = 2.0. But, the data 
do not show a scale invariance plot for different N which 
is expected from a scaling behavior. We have also im- 
plemented our numerical tool to calculate the exponent 
of energy gap at A = 0, 0.25, —0.25. Again, the plot 
of Ngap{N, h) versus Nh'^ shows a linear behavior for 
£ = 2.0. 

However, we should note that the horizontal axes pre- 
sented in Fig.f^l are limited to very small values of 
X = Nh^ < 0.0024. Thus, we are not allowed to read 
the scaling exponent of gap which exists in the thermo- 
dynamic limit {N — i- CX3 or X 3> 1). We have been limited 
to consider the maximum value of A^ = 24. Because for 
the present model (Eq.lQ)) the total S^ does not com- 
mute with the Hamiltonian and we should consider the 
full Hilbert space of 2^ in our computations. Moreover, 
to avoid the effect of level crossings, we should consider 
very small values oi h < 0.01. Therefore, the value of x 
can not be increased in this method. We will face up with 
the same problem even though the calculation is done 
by Density Matrix Renormalization Group (DMRG). In 
that case we can extend the calculation for larger sizes, 
A^ ~ 100 but the first level crossing happens for smaller 
value of h. The position of the first level crossing is pro- 
portional to 7V~^/^ which will be explained in the next 
sections. This is the level crossing between the higher 
excited states, in this case between the second and third 
ones. Thus, we have to find the scaling behavior from 
the small x data. 





FIG. 2: The value of gi{N) versus the chain length (N). (a) 
A = 0.5, (b) A = —0.5. Data for different transverse fields 
0.004 <h< 0.01 fall exactly on each other. 



Since g{N) is a function of only N we end up with 

e(f>s = 2. (21) 



IV. PERTURBATIVE APPROACH 

According to our computations where A^ < 24, the 
small X regime is equivalent to the very small h values. 
In this case, the energy gap of the finite size system is 
basically representing the perturbative behavior. Thus, 
Ea. (|15|l is rewritten in the following form 



gap{N, h)^^ + ^f{x) ^ ^ + g{N)h^ 



to a very good approximation. Because, the first order 
perturbation correction is zero in the transverse field and 
the leading nonzero term is /i^. In this way we get, 
g{N) — j^j^i . If the small x behavior of the scaling 



Nh'-' ■ 
function defined by: f{x) 



we find that, 



g{N) = B7v(^=-i)/i('^'^=-2). 



(20) 



The multiplication of both sides of Ea. H19|) by A^ leads to: 
Ngap{N, h) = B + Bx"^/" . This shows that in the small 
X regime, Ngap{N, h) is a linear function of x^/^. This is 
in agreement with our data in Fig.l^ where 4>s = 1 and 
according to Ea. (|21|) . the value of e is found to be e = 2. 
The function g{N) in Ea. (|19() is actually the coefficient 
of the first nonzero correction in the perturbation expan- 
sion for the energy gap of a finite chain. 



(19) gap{N,h) = gap{N,0) + gi{N)h^ + ■ ■ ■ + gm{N)h:' 



2m 



= §+9i{N)h^ 



9m{N)k 



2m 



(22) 



where m is an integer. The effect of higher order terms 
can be neglected for h < 0.01 to a very good approxi- 
mation. Now, let consider that the large N behavior of 
9i{N) is, 



N 



lim gi{N) ^aiN". 



(23) 



We find that, 



B 



gap{N,h)~-(l + b^N^+^h^), (24) 

where bi = ^ — const. We can write Ea. (|24(l in terms 
of the scahng variable {x = Nh'^), 



gap{N, h) 



B_ 

N 



l + A^^+i-i X 



(25) 



For large N limit, Eq.jSsJ should be independent of N. 
This impose the following relation, 

2 
a + l-- = 0. (26) 

£ 

This equation defines the relation between a and e, 

2 



a + 1 



(27) 



The above arguments propose to look for the large N 
behavior of gi(N). For this purpose, we have plotted in 
Fig-® the following expression versus N, 



9iiN) 



gapiN,h)-gap{N,0) 



(28) 



for fixed values of h (0.001 < h < 0.01) and A = 
0.5, —0.5. The results have been plotted for different 
sizes, N — 12, 14, ... , 24 to derive the a exponent de- 
fined in Ea. H23l) . In Fig.l^J-a) we have considered the 
case of A = 0.5 and found the best fit to our data for 
a = 0.88 ± 0.02. Therefore, e = 1.06 ± 0.02 which shows 
a very good agreement with Ea. lfTUI) . sft = y = 12/11 = 
1.09. Moreover, our data for different h values fall per- 
fectly on each other which show that our results for gi [N) 
is independent of h as we have expected. We have also 
plotted the results for A = —0.5 in Fig.|(21-b) and found 
a = 0.55 ± 0.02. Then, we obtain e = 1.29 ± 0.02 which 
can be compared with Eq.Q, Eft = M = 6/5 = 1.2. 
This shows a slight deviation which is the result of nu- 
merical computations and also the limitation on the size 
of system. Moreover, the magnitude of the energy gap 
is smaller for A — —0.5 than the case of A = 0.5 which 
implies less accuracy. However, the results for different 
h, give once more a unique gi{N) in agreement with its 
definition to be independent of h. 

We have extended our numerical computations to con- 
sider the other values of A. The results have been pre- 
sented in Table. I. We have listed a, the resulting e which 
is obtained from Eq. (|27l) and the result of field theoretical 
approach (sft) for different values of A. Our numerical 
results show very good agreement with the exponents de- 
rived in the field theoretical approacbiS. 



V. DISCUSSION AND SUMMARY 

We have studied the opening of excitation gap of the 
XXZ chain by breaking the rotational symmetry. We 



TABLE I: The a exponent defined in Eg. 11231 . the related gap 
exponent (e) and the corresponding value eft obtained by 
field theoretical approach for different anisotropy parameter 

(A). 



A 


a 


e 


epT 


0.70 


0.96 


1.02 


1.04 


0.50 


0.88 


1.06 


1.09 


0.25 


0.70 


1.17 


1.18 


0.0 


0.65 


1.21 


1.33 


-0.25 


0.60 


1.25 


1.26 


-0.50 


0.55 


1.29 


1.20 


-0.70 


0.59 


1.25 


1.14 



have implemented the modified Lanczos method to get 
the excited state energies at the same accuracy of the 
ground state one. The convergence of this method is fast 
for the low energy states but gets slower for higher energy 
ones. The reason is related to the condensation of states 
in a fixed interval of energy by going to higher states. 
Thus the method will be very sensitive to the initial pa- 
rameter Et (Eq.®). However, for the 10th lowest energy 
states we got fast convergence. We have been limited to 
consider the maximum N = 24, because for the present 
model (Eq.®) the total S'^ does not commute with the 
Hamiltonian. Thus, we should consider the full Hilbert 
space of 2^ in our computations. 

We have tried to find the scaling of energy gap in the 
presence of the transverse field by introducing the scal- 
ing function in Ea. H13() . According to this approach, the 
energy gap scales as gap ^ ft,^ where e defines x — Nh'^, 
the scaling variable. The right scaling exponent gives a 
linear behavior of Ngap(N, h) versus x for large x. To 
find the large x behavior we have faced with a serious 
problem. To see the scaling behavior we have to consider 
very small values of h to avoid the effect of level crossing 
between the excited state which defines the gap and the 
upper states. The position of the first level crossing is 
roughly proportional to the inverse of the chain length. 
Thus, we were not able to get the large x behavior and 
being far from the crossing point, at the same time. 

The limitation to very small h values states that our 
results are representing the perturbative ones. We have 
then led to get the scaling behavior by finding the di- 
vergence of gi{N) '^ N"' in the thermodynamic linrit 
{N -^ oo). The function gi{N) is the coefficient of the 
leading term in the perturbation expansion fEa. (|22() ). 
Based on the formulation presented in the previous sec- 
tion the gap exponent is related to the divergence of 
gi{N) by e = ^jrrj-- Our numerical results presented in 
Fig.® show that gi{N) is independent of h and its di- 
vergence versus N is defined by Ea. (|23|l . Moreover, the 
gap exponent which has been obtained by our numerics 
(listed in Table. I) is in very good agreement with the 
results obtained by the field theoretical approacbiSi. 



We might also find the gap exponent by finding the pre- 
ceise location of the first level crossing between the sec- 
ond and third excited states which we call hi. The slope 
of log{hi) versus log{N) is — . In the small h regime 
where Eq. (|22() is approximated up to the quadratic term, 
we can write two different expressions for the energy of 
the third and second excited states, namely 

E:i{N, h) - Eo{N, h) - gap^{N^ h^Q)+ g^i\N)h\ 
E2{N, h) - Eo{N, h) = gap2{N, h = Q)+ gf-\N)h\2Q) 

The subtraction of the two terms give the difference of 
E^ — E2 which is a function of the terms presented in 
the right sides of Eq. H29|l . The zero field terms are pro- 



portional to N~^ and the coefficients of the quadratic 
terms obey Ea. (|23|l . Thus, the scaling of hi in terms of 
the lattice size is like, hi ~ N^^^^ . However, the preceise 
determination of hi defines the accuracy of this approach. 
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